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ABSTRACT 


Steam  leakage  from  an  aircraft  carrier  catapult  is  sometimes  ingested  into  the 
aircraft  engines  upon  launch  which  may  induce  compressor  stall.  Investigation  of  this 
phenomenon  is  of  particular  interest  to  the  Navy  with  its  new  F35C,  the  aircraft  carrier 
variant  of  the  joint  strike  fighter.  The  single  engine  configuration  of  the  F-35C  makes 
this  aircraft  particularly  vulnerable  to  steam-induced  stall.  The  present  study  examined 
both  throttle-induced  stall  and  steam-induced  stall  in  a  compressor  at  90%  and  95%  speed 
through  the  use  of  9  Kulite  and  2  hot-film  pressure  transducers.  The  use  of  Fast  Fourier 
Transform  waterfall  plots  of  the  transient  data  before  and  during  stall  proved  invaluable 
in  determining  stall  precursors  as  well  as  the  mode  of  rotor  stall. 

In  addition,  a  new  computational  fluid  dynamic  model  was  designed  using  CFX-5 
software  to  represent  a  single  blade  passage  of  the  compressor  rotor,  in  order  to  predict 
compressor  performance.  The  computed  results  were  compared  to  experimental  results 
gathered  at  various  throttle  settings.  An  accurate  model  will  enable  researchers  to  predict 
compressor  performance  for  various  and  multiple  gases. 
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I.  INTRODUCTION 


As  the  U.S.  Navy  transitions  to  the  new  F-35C  Joint  Strike  Fighter  (JSF),  concern 
regarding  the  occurrence  of  steam-induced  stall  during  catapult  launches  has  increased. 
Research  conducted  at  Naval  Air  Engineering  Station  Lakehurst  with  an  F-18  has  shown 
the  susceptibility  of  jet  engines  to  stall  after  ingesting  steam.  The  F-35C,  with  a  single 
engine,  is  clearly  particularly  vulnerable  to  catastrophic  loss  should  the  engine  stall  as  a 
result  of  ingesting  steam  leaking  from  aging  aircraft  carrier  catapults. 

The  research  conducted  at  the  Turbopropulsion  Laboratory  at  the  Naval 
Postgraduate  School  is  presently  focused  on  the  effects  of  steam  ingestion  on  a  transonic 
compressor  rotor.  The  complete  rotor-stator  transonic  compressor  fan  stage  was 
specifically  designed  for  the  Naval  Postgraduate  School  by  Sanger  at  the  NASA  Glenn 
Research  Center  (Refs.  1  and  2),  to  be  evaluated  in  the  Transonic  Compressor  Rig  (TCR). 

The  performance  characteristics  of  the  compressor  in  both  the  fan-stage  as  well  as 
the  rotor-only  configurations  were  mapped  by  Gannon,  Hobson  and  Shreeve.  Utilizing 
this  data,  the  performance  characteristics  at  70%,  80%,  90%  and  100%  speed  prior  to  and 
during  stall  were  established  (Refs.  3-5).  Rodgers  reestablished  unsteady  pressure 
measurements  at  60%,  70%,  and  80%  speed  (Ref.  6).  Inlet  and  exit  surveys  at  70%,  80%, 
90%  and  100%  speed  were  conducted  with  a  three-hole  probe  by  Villescas,  who  then 
determined  the  spanwise  diffusion-factor  distribution  of  the  rotor  at  choke,  peak 
efficiency,  and  stall  (Ref.  7).  Brunner  repeated  the  inlet  surveys  with  a  5-hole  probe  and 
determined  the  pitch  angle  and  Mach  number  distributions  for  the  rotor  inlet  (Ref.  8). 
Payne,  using  both  hot-film  and  Kulite  pressure  transducers,  mapped  95%  design  speed  to 
near  stall  and  70%  design  speed  to  steam-induced  stall  (Ref.  9).  Levis  increased  the 
number  of  Kulite  pressure  transducers  to  nine,  mapped  90%  design  speed  up  to  steam- 
induced  stall,  and  re-established  the  70%  design  speed  stall  and  steam  stall  points  after 
the  replacement  of  the  inlet  honeycomb  flow  straightener  (Ref.  10). 

The  current  study  completed  the  feasible  perfonnance  mapping  of  the  rotor-only 
configuration  to  95%  of  design  speed.  By  replacing  the  previous  pneumatic  Scanivalve 
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system  with  three  Model  3217  Digital  Sensor  Arrays  (DSA)  and  rewriting  the  HP-VEE 
performance  program  to  interact  with  them,  pressure  readings  were  made  more 
efficiently,  which  resulted  in  less  wear  and  tear  on  the  TCR.  Nine  Kulite  pressure 
transducers  in  the  case  wall,  one  impact  Kulite  probe,  and  two  hot-film  transducers  were 
utilized  to  measure  rotor  characteristics  prior  to  and  during  steam-induced  stall.  The 
effect  of  steam  ingestion  upon  the  displacement  of  the  stall  line  at  90%  and  95%  speed, 
during  which  time  the  rotor  tip  relative  speeds  were  transonic,  was  measured. 

A  refined  computational  grid  (Ref.  11)  was  generated  of  a  single  rotor  passage  for 
the  numerical  investigation  of  the  steady-state  performance  of  the  rotor.  Favorable 
comparisons  with  the  experimental  data  at  90%  and  95%  speed  were  calculated. 
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II.  EXPERIMENTAL  FACILITY  AND  EQUIPMENT 


A.  TRANSONIC  COMPRESSOR  ROTOR 

For  the  most  recent  studies,  a  rotor-only  assembly  was  evaluated  in  order  to 
simplify  the  configuration  and  minimize  the  potential  damage  in  the  event  of  a 
catastrophic  failure  during  stall. 

The  Transonic  Compressor  Rig  (TCR),  as  shown  in  Figure  1,  was  driven  by  two 
opposed-rotor  air  turbine  stages,  supplied  by  a  12  stage  Allis-Chalmers  axial  compressor. 
The  rotor  had  22  blades  and  was  made  from  a  high  strength  aluminum  alloy  (7075-T6). 
The  design  specifications  for  the  Sanger  Stage  are  given  in  Table  1.  Air  was  drawn  into 
the  TCR  from  the  atmosphere  through  a  throttle  valve  and  a  honeycomb  flow 
straightener.  The  air  passed  down  a  five-meter  long  46  cm  diameter  pipe  connecting  the 
settling  chamber  to  the  test  compressor.  The  air  passed  through  a  nozzle,  used  to 
measure  the  mass  flow  rate,  and  was  exhausted  back  to  the  atmosphere  aft  of  the  test 
compressor. 


Figure  1 .  Transonic  Compressor  Rig,  with  Inlet  Piping  Removed 


3 


Table  1.  Sanger  Stage  Parameters 


Parameter 

Value 

Pressure  Ratio 

1.61 

Tip  Speed 

33.02  ms  (1300ft  s) 

Design  Speed 

27085  rpm 

Design  Mass  Flow 

7.75  kg'  s  (1 7.051b. ’s) 

Specific  Mass  Flow 

170.88  kg  s-nf  (351bm  s-fF) 

Specific  Head  Rise 

.246 

Tip  Inlet  Relative  Mach  Number 

1.28 

Aspect  Ratio 

1.2 

Hub  Tip  Radius  Ratio 

.51 

Rotor  Inlet  Ramp  Angle 

28.2 

Number  of  Rotor  Blades 

22 

Tip  Solidity  (Rotor) 

1.3 

Outside  Diameter 

0.2794m  (11  inches) 

Rotor  Diffusion  Factor  -  Tip 

.4 

Rotor  Diffusion  Factor  -  Hub 

.47 

B.  STEAM  INGESTION  SYSTEM 

In  order  to  generate  the  required  pressure  of  steam  for  this  experiment,  a  Sussman 
model  SVS600  steam  boiler,  capable  of  producing  saturated  steam  up  to  a  pressure  of 
1000  kPa  or  1.4  kg/sec  at  100  °C,  was  utilized  (Ref.  12).  The  transient  response  of  the 
steam  pressure  was  monitored  through  the  use  of  a  pressure  transducer  installed  in  the 
steam  pipe  (Figure  2).  Two  remotely-operated,  fast-acting  solenoid  valves  were  used  to 
release  the  steam  into  the  intake  plenum  as  well  as  to  initially  vent  the  pipe  to  preheat  and 
prevent  condensation.  Consistency  in  temperature  readings  from  the  inlet  thermocouple 
after  it  was  moved  from  the  6  o’clock  (Refs.  9  and  10)  to  the  8  o’clock  position 
confirmed  that,  after  ingestion  of  steam,  the  entire  compressor  inlet  saw  a  consistent 
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mixture  of  steam  and  air.  The  compressor  rig  and  steam  ingestion  system  are  shown 
schematically  in  Figure  2  and  a  photo  of  the  intake  plenum,  TCR,  inlet  piping,  and  steam 
generator  are  shown  in  Figure  3. 


Kulite  pressure 
tranduscers 
(3  circumferentially 
followed  by  6  axially 
in  the  casewall) 


Figure  2.  Transonic  Compressor  Rig  with  Steam  Ingestion  System 
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Figure  3.  Intake  Plenum  and  Steam  Ingestion  System 


6 


III.  INSTRUMENTATION 


A.  PROBES 

1.  Hot-Film 

Inlet-flow  conditions  were  characterized  using  hot-film  anemometry 
measurements.  Measurements  were  taken  using  two  20  micron  TSI  model  1212-20  hot- 
film  probes  (Figure  4,  Ref.  12)  with  serial  numbers  70618113  and  70618114.  Data  were 
acquired  using  an  IFA  100  Intelligent  Flow  Analyzer,  and  TSI  ThennalPro  software 
version  4.53,  running  on  a  PC. 


7.9  mm  (..112) 
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3.2  mm  (.125)  Dia.  4.G  mm  (.IQ.)  Did 


Figure  4.  TSI  Model  1212-20  Hot-Film  Probe  (From  Ref.  12) 


2.  Kulite  Pressure  Transducer 

A  model  XCQ-080-25  Kulite  miniature  silicon  pressure  transducer  was  used  to 
obtain  time-resolved  pressure  data.  The  probe  was  a  miniature,  semiconductor,  strain- 
gauge  transducer  which  incorporated  a  fully  active  four-arm  Wheatstone  bridge 
dielectrically  isolated  silicon-on-silicon  diaphragm.  A  diagram  of  the  probe  is  given  in 
Figure  5  and  the  specifications  are  given  in  Table  2. 


Figure  5.  Kulite  XCQ-080  Series  Transducer  (From  Ref.  6) 
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Table  2.  XCQ-080-25  Factory  Specifications  (From  Ref.  5) 


Input 

Pressure  Range 

25  psi 

Over  Pressure 

50  psi 

Burst 

75  psi 

Rated  Electrical  Excitation 

10VDC/AC 

Maximum  Electrical  Excitation 

15VDC/AC 

Input  Impedance 

800  Ohms 

Output 

Output  Impedance 

1000  Ohms 

Full  Scale  Output 

100  mV 

Residual  Unbalance 

+-3%  FSO 

Non-Linearity  and  Hysterisis 

0.1  %  FS  BFSL 

Hysteresis 

0.1% 

Repeatability 

0.1% 

Resolution 

Infinite 

Natural  Frequency 

300kHz 

Perpendicular  Accel  Sensitivity 

0.0003%  FS/g 

Transverse  Accel  Sensitivity 

0.00004%  FS/g 

Insulation  Resistance 

100  Megohm 

Environmental 

Operating  Temp  Range 

-65  to  250  deg  F 

Compensated  Temp  Range 

80  to  180  deg  F 

Thermal  Zero  Shift 

+- 1  %  FS/100  F 

Thermal  Sensitivity  Shift 

+- 1  %  FS/100  F 

3.  Fine-Wire  Thermocouple 

An  OMEGA  fine-gage  thermocouple  was  used  to  obtain  fast,  accurate 
temperature  measurements  of  the  steam  in  the  inlet  to  the  TCR.  The  thermocouple 
utilized  J-Type  wires  of  diameter  0.025  mm  and  was  capable  of  obtaining  accurate 
temperature  measurements  in  both  air  and  water.  A  diagram  of  the  thermocouple  is  given 
in  Figure  6  and  the  specifications  are  given  in  Table  3. 
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Figure  6.  Fine  Gage  Microtemp  Thermocouple 

Table  3.  Thennocouple  Response  Time 


Wire  Size 
mm  (inches) 

Still  Air 
427C/38C 
800F/100F 

60  ft/sec  Air 
427C/38C 
800F/100F 

Still  H20 
93C/38C 
200F/100F 

0.025  (0.001) 

0.05  sec 

0.004  sec 

0.002  sec 

B.  INSTALLATION  OF  SENSORS 
1.  Hot-Film 

Two  hot-film  probes  were  installed,  as  shown  in  Figure  7,  through  the  case  wall 
1 1.5cm  upstream  of  the  leading  edge  of  the  rotor,  and  extended  4cm  into  the  flow.  The 
sensors  were  positioned  at  the  2  and  10  o’clock  positions  and  faced  toward  the  inlet  flow. 
Adding  a  second  probe  allowed  for  redundancy,  in  case  of  failure,  and  for  comparison  of 
results. 
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Figure  7.  Probe  Location  in  TCR  Inlet 

2.  Kulite  Pressure  Transducer 

Nine  Kulite  pressure  transducers  were  placed  in  the  case  wall  at  various  locations 
axially  and  peripherally  to  examine  the  effect  of  stall  on  the  case  wall  pressure,  as  shown 
in  Figures  8  and  9. 
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One  impact  Kulite  probe  was  installed  at  the  6  o’clock  position.  The  probe  was 
directed  upstream  and  had  a  vertical  row  of  holes  drilled  on  its  face  instead  of  the 
standard  screen  which  is  shown  in  Figure  5. 

C.  DATA  ACQUISITION 

1.  Hot-Film  and  Kulite  Data  Acquisition 

Hot-film  data  were  acquired  as  described  by  Payne  (Ref.  9).  Kulite  Data  were 
acquired  as  described  in  detail  by  Levis  (Ref.  10).  Figure  10  is  the  combined  Kulite  and 
hot-film  data  acquisition  system. 
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Figure  10.  Transient  Hot-Film  and  Kulite  DAC  Systems 

2.  Digital  Sensor  Array 

The  Digital  Sensor  Array  is  a  stand-alone  temperature  compensated  electronic 
pressure  scanner  capable  of  up  to  16  individual  or  8  differential  pressure  readings  (Ref. 
13).  Each  pressure  sensor  is  calibrated  to  work  within  a  temperature  range  of  0  to  59°C. 
Each  unit  is  assigned  a  specific  IP  Address  enabling  it  to  receive  commands  and  return 
data  via  the  Ethernet.  The  unit  IP  Address  had  to  be  renamed  in  order  to  make  it 
compatible  with  the  laboratory  local  area  network  (LAN)  (see  Appendix  A). 

Data  were  transferred  directly  from  the  DSA  3217  pressure  bricks  to  an  EXCEL 
spreadsheet  via  the  Ethernet  connection  and  the  HPVEE  computer  program  TCR.vee, 
which  converted  the  binary  data  to  engineering  units  (Pa)  (Refs.  13  and  14).  Pressure 
readings  could  also  be  acquired  individually  from  the  bricks  using  the  DSALink  3 
program,  Figure  11,  by  pressing  the  start  button  (traffic  sign)  and  the  stop  button  (stop 
sign). 
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IV.  EXPERIMENTAL  PROCEDURE 


A.  HOT-FILM  CALIBRATION 

The  hot-film  probes  were  individually  calibrated  after  being  positioned  inside  the 
compressor  inlet,  using  the  calibration  procedure  detailed  in  Appendix  B.  First  a  shorting 
probe  was  inserted  and  the  cable  resistance  recorded  with  the  IF  A  100.  The  shorting 
probe  was  removed,  replaced  by  a  hot-film  probe,  and  the  new  resistance  read.  The  IFA 
100  was  adjusted  for  the  probe’s  operating  resistance  and  the  appropriate  bridge 
compensation.  The  compressor  was  started  and  its  speed  increased  to  30%  of  maximum 
speed.  Once  speed  stabilized,  the  differential  pressure,  the  difference  between  static  and 
stagnation  pressure,  was  calculated  and  entered  into  the  ThermalPro  calibration  software, 
and  a  calibration  point  was  recorded.  This  step  was  repeated  in  increasing  increments  of 
10%  to  100%  speed.  The  software  applied  curve  fitting  and  calculated  King’s  Law 
coefficients  which  were  recorded  and  saved  for  future  data  analysis.  During  this  process, 
a  constant  density  was  assumed  which  led  to  a  maximum  possible  error  of  12%. 

B.  DIGITAL  SENSOR  ARRAY  SETUP  AND  OPERATION 

The  DSA  units  were  factory  calibrated,  so  no  calibration  was  required.  The  units 
were  energized  and  at  least  30  minutes  was  allowed  for  warm  up  prior  to  use.  The 
DSALink  3  software  was  used  to  connect  to  each  unit  to  check  its  status  and  conduct  a 
zero  calibration  prior  to  pressures  being  applied.  The  ScanPressures.vee  program  was 
run  to  complete  DSA  setup,  changing  its  data  output  to  binary  and  pressure  readings  to 
Pascals.  The  TCR.vee  program  (Figure  12)  was  used  during  compressor  operation  to 
measure  steady-state  pressures  and  temperatures  throughout  the  TCR  and  calculate  the 
performance  of  the  compressor  rotor. 
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Figure  12.  TCR.Vee  Program 


C.  COMPRESSOR  OPERATION 

The  standard  procedure  for  operating  the  TCR  was  to  increase  to  the  appropriate 
speed  with  the  throttle  fully  open.  While  maintaining  constant  rotor  speed,  the 
compressor  throttle  position  was  changed  so  measurements  could  be  taken  at  different 
mass  flow  rates.  By  measuring  mass  flow  rate,  inlet  and  exit  total  temperatures  and 
pressures  the  total-to-total  pressure  ratio  and  isentropic  efficiency  were  calculated  and  the 
compressor  performance  was  mapped  for  various  speeds  from  open  throttle  to  throttle- 
induced  stall.  Measurement  and  calibration  procedures  were  described  in  more  detail  by 
Gannon  et  al.  (Ref.  3).  Normal  speeds  used  are  shown  in  Table  4. 


Table  4.  TCR  Speed  Settings 


Percent  Speed 

RPM 

100 

27085 

95 

25730 

90 

24376.5 

70 

18959.5 
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D.  HOT-FILM  MEASUREMENTS 

Data  were  taken  at  90  and  95%  speed  in  order  to  calculate  the  compressor  inlet 
average  velocity  and  turbulence  intensity.  Hot-film  measurements  were  also  used  to 
study  the  airflow  prior  to  and  during  stall  to  investigate  possible  precursors  to  stall. 
Measurements  were  taken  at  open  throttle,  peak  efficiency,  near  stall,  throttle-induced 
stall,  and  steam-induced  stall.  The  steps  in  using  the  IFA  100  to  collect  data  for  two  hot- 
film  probes  are  detailed  in  Appendix  C.  The  ThermalPro  software  was  set  to  read  both 
probes  at  a  rate  of  100,000  Hertz.  For  open  throttle,  peak  efficiency,  and  near  stall,  the 
sample  size  was  set  to  8  kpts/ch,  which  lasted  0.0819  seconds.  The  sample  size  was 
increased  to  5 12  kpts/ch  for  the  stall  runs  to  ensure  that  both  the  precursor  and  stall  event 
were  captured  in  a  single  file  which  lasted  5.2429  seconds. 

E.  STALL  RUNS 

In  order  to  initiate  a  throttle-induced  stall,  the  compressor  speed  was  increased  to 
90  or  95%  speed.  The  inlet  throttle  was  closed  incrementally,  with  readings  taken  at  each 
setting,  until  compressor  stall  occurred.  For  both  speeds,  stall  occurred  as  the  throttle 
was  closed  from  the  7.2  position  (72%  closed). 

For  steam-induced  stall,  the  compressor  was  throttled  to  near  stall  conditions,  at 
which  time  steady-state  measurements  were  taken.  The  boiler  isolation  valve  and  steam 
vent  solenoid  valve  were  both  opened  to  heat  the  steam  piping  in  order  to  prevent 
condensation.  After  the  piping  was  fully  heated,  the  vent  valve  was  closed  and  the  data 
trace  from  the  pressure  transducer  and  thermocouple  was  started,  (Figure  13).  When  the 
pressure  in  the  pipe  reached  a  target  pressure,  the  isolation  valve  was  closed.  A 
countdown  to  steam  release  commenced,  at  which  time  the  Kulite  and  hot-film  data 
acquisition  were  initiated.  At  the  end  of  the  countdown,  the  main  steam  valve  (fast¬ 
acting  solenoid  valve)  was  opened,  allowing  the  steam  to  dump  into  the  inlet  plenum  of 
the  compressor.  The  suction  at  the  inlet  drew  the  steam  down  the  piping  and  into  the 
compressor.  After  several  seconds,  all  data  acquisition  was  stopped.  The  whole 
procedure  was  repeated  at  reduced  throttle  settings  until  stall  occurred. 
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Figure  13.  Strip  Chart  Program 
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V.  PROCEDURE  FOR  COMPUTER  SIMULATION 


A.  ICEM-CFD 

ICEM-CFD  was  a  CFX-5  compatible  Computer  Aided  Design  (CAD)  program 
that  allowed  a  user  to  create  or  import  the  geometry  of  their  choice.  In  this  case,  the 
imported  geometry  (Ref.  15)  was  of  a  single  blade  and  its  passage  from  inlet  to  outlet  of 
the  compressor.  A  block  was  created  to  encompass  the  entire  blade  passage,  its  nodes 
made  periodic  (circumferential  angle  16.3636°),  and  its  edges  associated  to  the  passage 
edges.  An  H-grid  was  created  around  the  blade  and  its  edges  and  vertices  associated  to 
those  of  the  blade.  Extra  vertices  were  added  to  fit  the  H-grid  to  the  curvature  of  the 
blade.  After  selecting  the  block  representing  the  blade  only,  an  O-grid  was  formed  to 
enable  the  user  to  refine  the  mesh  near  the  blade.  The  block  representing  the  blade  was 
then  deleted  in  order  to  create  a  mesh  for  only  space  in  which  air  could  flow.  After 
creating  a  mesh  over  the  entire  grid  with  a  maximum  element  size  of  0. 125,  the  mesh  was 
refined  by  increasing  the  number  of  nodes  to  gamer  greater  data  regarding  the  tips  and 
surface  of  the  blade  than  the  inlet  or  outlet.  Figures  14,  15,  and  16  show  different  views 
of  the  final  mesh  created  and  Table  5  gives  the  final  number  of  elements  used  to  generate 
the  grid. 


Figure  15.  Hub  Assembly  Mesh  Created  in  ICEM-CFD 


Figure  16.  Shroud  Mesh  Created  in  ICEM-CFD 
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Table  5.  Number  of  Elements  within  the  Mesh 


Part 

Elements 

Blade 

3280 

Curves 

1035 

Exit 

1000 

Hub 

4293 

Inlet 

1000 

Passage  1 

5560 

Passage  2 

5560 

Points 

14 

Solid 

171720 

Tip 

4293 

Total  Elements 

197757 

Total  Nodes 

184418 

B.  CFX-5 

Previous  studies  capable  of  creating  fairly  accurate  computational  models  of  the 
rotor-only  configuration  were  conducted  with  the  SWIFT  program  using  a  very  refined 
C-grid  (Ref.  5).  Unfortunately,  the  SWIFT  program  was  not  capable  of  incorporating  a 
steam/air  mixture.  The  new  challenge  has  been  to  create  a  refined  grid  in  the  CFX-5 
program  which  had  a  function  to  input  a  variety  or  mixture  of  gases  and  multiphase 
flows. 


CFX-5  is  comprised  of  three  main  functions:  CFX-Pre,  CFX-Solver,  and  CFX- 
Post.  In  CFX-Pre,  the  “turbo”  mode  was  specified  to  ensure  the  program  modeled 
turbomachinery  flows.  To  most  accurately  model  the  real  compressor,  a  k-s  turbulence 
model  was  used  with  the  values  of  k  and  s  calculated  from  the  experimentally  calculated 
values  of  average  velocity  and  turbulence  intensity.  Varied  average  exit  static  pressures 
were  input  as  boundary  conditions  to  model  the  throttling  of  the  compressor.  Figure  17 
shows  the  model  after  all  initial  values  were  set.  Further  details  regarding  CFX-Pre  setup 
may  be  found  in  Bochette  (Ref.  1 1).  CFX-Solver  was  set  to  conduct  1000  iterations  in  an 
attempt  to  converge  the  residual  target  to  lxlO"6.  If  the  solution  was  unable  to  converge 
in  the  first  1000  iterations,  the  CFX-Solver  program  allowed  the  user  to  continue  solving 
from  a  previous  file  run.  In  CFX-Post,  the  results  could  be  reviewed  visually  with  the 
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turbo  plotter  which  created  color  images  of  the  blade  passage  and  any  variable  the  user 
desired.  The  program’s  function  calculator  was  used  to  calculate  the  mass  flow,  mass- 
averaged  pressure  ratio,  and  temperature  ratio  to  facilitate  efficiency  calculations. 


Figure  17.  CFX-Pre  Model  of  a  Single  Blade  Passage  of  the  TCR 
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VI.  RESULTS  AND  DISCUSSION  OF  COMPRESSOR  MAPPING 


A.  COMPRESSOR  CHARACTERISTICS 

Average  velocity  and  turbulence  intensity  values  were  calculated  using  MATLAB 
code  “hw.m”  shown  in  Appendix  D  (Ref.  16).  These  values  versus  mass  flow  were  then 
plotted  in  an  EXCEL  spreadsheet,  as  seen  in  Figures  18  through  21. 


90%  Average  Velocity 


Figure  18.  90%  Speed  Average  Velocity  vs.  Mass  Flow  Rate 

90%  Speed  Turbulence  Intensity 
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Figure  19.  90%  Speed  Turbulence  Intensity  vs.  Mass  Flow  Rate 
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95%  Speed  Average  Velocity 
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Figure  20.  95%  Speed  Average  Velocity  vs.  Mass  Flow  Rate 

95%  Speed  Turbulence  Intensity 
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Figure  21.  95%  Speed  Turbulence  Intensity  vs.  Mass  Flow  Rate 


The  average  inlet  velocity  increased  quadratically  with  throttle  setting  for  both  90 
and  95%  speed.  At  90%  speed,  the  average  turbulence  intensity  was  2.3%.  At  95%,  the 
average  turbulence  intensity  was  2.2%.  As  may  have  been  expected,  this  compressor 
characteristic  appeared  to  be  independent  of  compressor  speed  and  throttle  setting. 
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B.  COMPRESSOR  SPEED  PERFORMANCE  DATA 

Continuing  on  compressor  testing  conducted  by  Gannon  et  ah,  Payne,  and  Levis 
(Refs.  3,  9,  and  10  respectively),  performance  data  from  open  throttle  to  stall  were 
obtained.  The  results  can  be  seen  in  Figures  22  and  23. 
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Figure  22.  Compressor  Pressure  Ratio 


Figure  23.  Compressor  Efficiency  (Fourth  Order  Polynomial  Trendlines) 
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As  can  be  seen,  the  stall  line  was  displaced  to  the  left  due  to  steam  ingestion  of 
0.12  kg/sec.  For  a  particular  speed  line,  the  mass  flow  rate  based  stall  margin  is  the  mass 
flow  rate  at  peak  efficiency  minus  that  at  stall  divided  by  the  peak  efficiency  mass  flow 
rate.  At  90%  speed,  mass  flow  rate  based  stall  margin  was  reduced  from  0. 17  to  0. 16  due 
to  a  2%  steam  mass  flow  rate.  At  95%  speed,  mass  flow  rate  based  stall  margin  was 
reduced  from  0.14  to  0.12. 


Steam  ingestion  resulted  in  an  increase  in  inlet  temperature,  which  caused  fluid 
density  to  decrease  and  the  rotor  to  overspeed.  This  effect  is  displayed  in  Figures  24  and 
25,  which  shows  the  drop  in  pressure  in  the  steam  line  and  the  subsequent  increase  in 
inlet  temperature  due  to  steam  ingestion  during  a  stall  event. 
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Figure  24.  Steam  Effects  at  90%  Speed 
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Figure  25.  Steam  Effects  at  95%  Speed 

C.  PEAK  EFFICIENCY  AND  NEAR  STALL  STEADY  STATE 

Hot-film  and  Kulite  measurements  were  made  at  peak  efficiency  and  near  stall. 
Contour  plots  and  auto-correlation  plots  were  developed  from  the  measurements  to 
develop  a  better  understanding  of  the  compressor  characteristics  and  inlet  flow  behavior 
at  different  mass  flow  rates. 

An  auto-correlation  plot  was  created  with  experimental  hot-film  data  and 
MATLAB  code  fftauto.m,  Appendix  D,  and  used  to  estimate  the  value  of  the  integral 
time  scale,  A  t.  The  area  under  the  curve  is  an  indication  of  the  length  scale  of  the  inlet 
turbulence.  Figure  26  illustrates  an  example  auto-correlation  of  data  at  90%  speed. 
Below  is  a  sample  calculation. 
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Figure  26.  Auto-Correlation  of  Hot-Film  Data  at  90%  Speed  Near  Stall 


At  *  0.002 
Ti  =  0.020333 
Uavg  =  62.9 

u’2  =  (Ti  *  Uavg)2  =  1.6357 
k  =  1.5*(u’2)  =  2.45355 
L  =  At  *  Uavg  =  0.1258 
8  =  (C°-75*k1-5)/(K  *  L)=  12.55 


(From  above  graph) 

(Turbulence  intensity  from  hwl.m) 
(Average  velocity  from  hwl  .m) 

(Turbulence  length) 

(C  =  0.09,  k  =  0.4) 


Contour  waterfall  plots  were  created  for  Kulite  #1,  the  impact  Kulite  and  hotwire 
#1  for  the  cases  of  peak  efficiency  and  near  stall  at  90%  speed,  Figures  27  through  32. 
The  process  was  repeated  for  95%  speed,  Appendix  E. 
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Figure  27.  Kulite  1  at  90%  Speed  Peak  Efficiency 


Figure  28.  Impact  Kulite  at  90%  Speed  Peak  Efficiency 
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The  hot-film  probes  did  not  detect  the  once-per-rev  signal  because  they  were  too 
far  upstream  of  the  rotor,  however  they  did  respond  to  the  precursor  signal  which  was  a 
flow  phenomenon.  The  Kulites  being  over  the  rotor  did  detect  a  once-per-rev  signal. 
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Figure  29.  FIot-Film  1  at  90%  Speed  Peak  Efficiency 
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Figure  30.  Kulite  1  at  90%  Speed  Near  Stall 
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Figure  3 1 .  Impact  Kulite  at  90%  Speed  Near  Stall 
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Figure  32.  Hot-Film  1  at  90%  Speed  Near  Stall 
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A  stall  precursor  was  only  evident  in  the  case-wall  Kulite  data  (Figure  30).  This 
was  not  consistent  with  data  taken  by  Payne  (Ref.  9)  who  detected  a  stall  precursor  with 


his  hot-film  probe.  An  explanation  for  this  discrepancy  could  be  that  the  hot-film  probe 


measurements  in  the  present  study  were  contaminated  due  to  probe  vibrations. 

D.  MEASUREMENTS  OF  90%  SPEED  STEAM-INDUCED  STALL 

Experimental  data  were  acquired  through  the  use  of  hot-film  probes  during  the 
steam-induced  stall  of  the  compressor  at  90%  speed.  A  power  spectrum  was  created 
using  a  Fast  Fourier  Transform  in  the  MATLAB  program  fftauto.m  (Appendix  D).  The 
result  gained  from  the  raw  data  (Figure  33)  can  be  seen  in  Figure  34.  The  FFT  displayed 
the  strong  effect  the  shock  waves  propagated  upstream  from  the  rotor  had  on  the  hot-film, 
causing  geometry  effects  like  once  per  revolution  frequency  (406  Hertz)  to  be  blocked, 
but  not  the  precursor  frequency  (approximately  250  Hertz). 
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Figure  33.  Hot-Film  Raw  Data  at  90%  Speed  During  a  Steam-Induced  Stall 
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Figure  34.  Hot-Film  Fast  Fourier  Transform  at  90%  Speed  During  a  Steam-Induced 

Stall 

The  contour  plot  of  the  hot-film,  Figure  35,  displayed  similar  effects,  making  the 
once  per  revolution  and  precursor  frequencies  difficult  to  distinguish.  The  Kulite 
pressure  transducers,  placed  in  the  case  wall  and  closer  to  the  rotor,  displayed  these 
frequencies  with  much  greater  distinction  as  seen  in  Figure  36.  This  plot  also  reinforced 
the  estimated  precursor  frequency  of  250  Hertz. 
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Figure  35.  Contour  Plot  of  Hot-Film  1  as  Compressor  Moves  into  a  90%  Speed  Steam- 

Induced  Stall 
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Figure  36. 
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Contour  Plot  of  Kulite  1  as  Compressor  Moves  into  a  90%  Speed  Steam- 

Induced  Stall 


E.  MEASUREMENTS  OF  95%  SPEED  STEAM-INDUCED  STALL 

The  procedure  described  in  part  D  was  repeated  in  the  TCR  for  95%  speed  steam- 
induced  stall  with  the  results  shown  in  Figures  37  to  40  and  additional  plots  in  Appendix 

F.  The  precursor  frequency  (approximately  300  Flertz)  and  its  harmonic  dominated  the 
FFT,  making  the  spike  expected  for  the  once  per  revolution  frequency  (429  Flertz) 
indiscernible. 


34 


Figure  37.  Hot-Film  Raw  Data  at  95%  Speed  During  a  Steam-Induced  Stall 


Figure  38.  Hot-Film  Fast  Fourier  Transform  at  95%  Speed  During  a  Steam-Induced 

Stall 
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Figure  39.  Contour  Plot  of  Hot-Film  1  as  Compressor  Moves  into  a  95%  Speed  Steam- 

Induced  Stall 
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Figure  40.  Contour  Plot  of  Kulite  1  as  Compressor  Moves  Into  a  95%  Speed  Steam- 

Induced  Stall 
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VII.  RESULTS  AND  DISCUSSION  OF  COMPUTER  SIMULATION 


A.  90%  SPEED  SIMULATION 

Using  the  average  velocity  and  k  and  s  values  calculated  from  experimental  data 
and  the  CFX-5  computer  software,  the  computer  simulated  rotor  performance  was 
mapped  from  open  throttle  to  near  stall  conditions  at  90%  speed.  Figure  41  displays  the 
pressure  on  the  blade  and  hub  and  the  Mach  number  distribution  at  95%  span.  Figure  42 
is  a  plot  of  the  static  pressure  over  the  surface  of  the  blade  at  its  mid-span  at  near  stall 
conditions. 


Mach  Number 

(Contour  1) 

II .  297etOOO 
1 . 197e+000 
t ,  097e+000 
9  975e- 00 1 
8  977e-001 
7 . 980e - 00 1 
6.982e-001 
5 . 985e - 00 1 
4  987e-001 
3  990e-001 
2 . 992e -00 1 
1 . 995e - 00 1 
I  9 . 975e-002 
*  1. 0006-015 


Pressure 

(Blade) 
3.522e+004  m 


-3.331e+004 


Figure  41 .  Computer  Simulation  Model  of  TCR  at  90%  Speed  Near  Stall 
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Figure  42.  Blade  Tip  Pressure  Distribution  at  90%  Speed  Near  Stall 


B.  95%  SPEED  SIMULATION 

The  process  described  in  part  B  was  repeated  for  95%  speed.  The  results  are 
displayed  in  Figures  43  and  44. 
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Figure  43.  Computer  Simulation  Model  of  TCR  at  95%  Speed  Near  Stall 
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Figure  44.  Blade  Tip  Pressure  Distribution  at  95%  Speed  Near  Stall 


C.  COMPARISON  OF  COMPUTED  AND  MEASURED  PERFORMANCE 

The  computer  simulated  values  of  mass  flow  rate,  total  pressure,  and  total 
temperature  were  used  to  calculate  pressure  ratio  and  isentropic  efficiency.  The 
spreadsheet  results  are  presented  in  Appendix  G.  Figure  45  is  a  graphical  comparison  of 
computer  simulated  and  experimental  pressure  ratio  results.  The  simulation  data  follows 
the  general  curvature  of  the  experimental  data.  However,  the  simulation  underestimated 
the  pressure  ratio  by  2.5  to  4.5%  as  it  approached  near  stall  conditions,  and  at  open 
throttle,  was  even  less  accurate.  The  simulation  also  reached  stall  at  a  higher  mass  flow 
than  was  found  experimentally. 


40 


Figure  45.  Experimental  and  Computer  Simulated  Pressure  Ratio  vs.  Mass  Flow  for  90 

and  95%  Speed  TCR 


Figure  46  shows  the  experimental  and  computer  simulated  efficiency  values.  As 
is  often  the  case  with  current  computer  models,  efficiency  was  overestimated.  In  this 
case,  peak  efficiency  was  overestimated  by  approximately  3.3%,  while  the  mass  flow  at 
which  the  peak  occurred  was  underestimated  by  approximately  6.5%. 
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Figure  46.  Experimental  and  Computer  Simulated  Efficiency  vs.  Mass  Flow  for  90  and 

95%  Speed  TCR 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  Compressor  Performance 

The  present  experiment  provided  evidence  of  the  detrimental  effect  steam 
ingestion  had  on  the  stall  line  of  the  compressor.  Stall  occurred  at  higher  mass  flow 
values  when  steam  was  ingested.  The  presence  of  steam  decreased  the  density  of  the 
fluid  flowing  through  the  compressor,  causing  the  compressor  to  overspeed,  then  stall. 
At  90%  speed  the  stall  margin  was  reduced  from  0.17  to  0.16  due  to  steam  ingestion  and 
at  95%  speed  it  was  reduced  from  0.14  to  0.12.  Approximately  2%  of  total  mass  flow 
rate  of  steam  was  ingested  by  the  rotor  for  both  these  speeds.  The  present  investigation 
was  the  first  transonic  tip  speed  steam-induced  stall  tests  on  the  TCR. 

Although  the  precursors  to  stall  were  as  distinctly  measured  in  the  present  study 
as  those  measured  previously  at  90%  and  95%  speed,  the  addition  of  multiple  Kulite 
probes  and  back-up  hot-film  probes  proved  invaluable  in  detennining  stall  precursors  at 
transonic  rotor  tip  speeds. 

2.  Computer  Simulation 

The  CFX-5  software  proved  to  be  a  relatively  accurate  modeling  system  for  the 
transonic  compressor.  It  will  prove  useful  for  modeling  the  effects  of  steam  ingestion  in 
a  rotor-only  compressor  as  well  as  a  single  or  multi-stage  compressor  at  near  stall 
conditions.  However  the  present  model  is  not  accurate  enough  to  use  for  design 
purposes. 

B.  RECOMMENDATIONS 

The  next  step  is  to  examine  the  effects  of  steam  ingestion  on  a  full  single  stage 
compressor  at  near  stall  conditions.  Steam-induced  surge  could  prove  catastrophic  so 
care  needs  to  be  taken  when  undertaking  such  tests. 
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The  computer  simulation  may  be  improved  by  further  refining  the  grid.  Modeling 
of  the  entire  compressor  instead  of  a  single  blade  passage  would  also  provide  a  more 
accurate  simulation  of  the  compressor.  With  experimental  steam  data  now  available  and 
the  CFX-5  software  capable  of  modeling  steam  as  the  fluid  medium,  a  model  of  the 
compressor  during  steam  ingestion  should  be  attempted. 
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APPENDIX  A.  RENAME  DSA  3217  IP  ADDRESS 


Step  1:  De-energize  DSA.  Connect  trigger/serial  test  cable  from  DSA  to  I/O  port  on 
computer. 

Step  2:  Create  a  HyperTenninal 

-On  computer  desktop  click  START— ♦PROGRAMS— ^ACCESSORIES 
—♦COMMUNICATIONS— ♦HYPERTERMINAL 
-Do  not  make  HyperTenninal  the  default.  Click  “NO” 

-Enter  a  name  and  pick  symbol 

-Connect  using  COM  1  -  or  whichever  I/O  the  serial  cable  is  connected  to 
-Set  Parameters  — ♦  9600  BAUD,  8  data  bits,  No  parity,  1  stop  bit 

Step  3:  Energize  DSA.  System  information  will  display. 

Step  4:  “Press  any  key  to  stop  auto-boot...”  displays.  Hit  any  key  within  3  seconds. 

Step  5:  When  the  auto-boot  process  is  stopped,  type  “p”  <enter>  to  get  existing 
parameters. 

Step  6:  Type  “c”  <enter>  to  change  parameters.  The  boot-loader  prompts  for  each 
parameter.  If  correct  press  <enter>. 

Step  7:  When  it  gets  to  “inet  on  ethernet,”  type  “131.120.20.###  (the  last  3  numbers  will 
be  between  250  and  300). 

Step  8:  After  all  changes  have  been  made,  verify  new  settings  by  typing  “p”  <enter>. 

Step  9:  Restart  the  operating  system  with  the  new  settings  by  typing  “@”  <enter>. 

Step  10:  If  DSA  boots  correctly,  de-energize  and  disconnect  the  serial  test  cable. 
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APPENDIX  B.  HOT-FILM  CALIBRATION  PROCEDURE 


IFA  100  Setup 

1 .  Record  current  temperature  (Ta)  and  pressure  (Pa). 

2.  Set  IFA  for  appropriate  channel  [CHANNEL] [#] [ENTER]. 

3.  Insert  shorting  probe  and  measure  cable  resistance  (Re). 

a.  Press  [RES  MEAS] 

b.  Zero  with  [OPERATE  RES]  control  knob  of  correct  transducer 

c.  Press  [RES  MEAS] 

d.  Record  Re  then  press  [ENTER] 

4.  Replace  the  shorting  probe  with  the  hot  wire  probe. 

5.  Repeat  step  #2  a  through  c  only.  DO  NOT  PRESS  ENTER.  Record  Ro. 

6.  Press  [OPERATE  RES]  and  adjust  to  hot  wire  operating  resistance  (on 
box  cover).  Record  Rop. 

7.  Press  [BRIDGE  COMP]  and  set  with  [BRIDGE]  control  knob  according 
to  following  table: 


Sensor  Type 

Standard  1 

T1.5 

35 

-P2 

83 

-PI  1.5 

115 

-PI5 

120 

Metal  Clad 

775 

-10A 

70 

-10 

70 

-20 

115 

-60 

250 

8.  Press  [RUN]. 

9.  Start  wind  tunnel.  Ensure  reference  pressure  is  set  to  10  inches  Hg. 
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10.  Turn  [CABLE  COMP]  knob  fully  counterclockwise.  Turn  knob  until 
OSC  light  comes  on  then  turn  clockwise  until  it  goes  off  plus  an  extra  half 
turn. 

11.  Create  a  new  folder  for  day.  Move  copy  of  calibration  file  70536123  into 
folder. 

12.  On  IFA  main  screen  select  [CALIBRATION]  then  [PROBE  DATA] 

a.  Select  [OPEN  CAL  FILE] 

b.  Select  file  70536123 

c.  Settings  as  follows: 


A/D  Chan 

Hotwire  channel  # 

IFA  Chan 

Hotwire  channel  # 

Probe  Type 

S 

Cable 

Re 

Opr  Res 

Rop 

Offset 

3 

Gain 

4 

Cal  Method 

Acquire  E  &  Type 
dP  Manual 

13.  Select  [CALIBRATE]. 

14.  Enter  current  conditions  on  the  left  hand  side  of  the  screen  remembering  to 
select  appropriate  dP  unit.  The  dP  is  the  difference  between  static  and 
stagnation  pressure  (Ports  6  +  48)  and  is  used  by  the  software  to  calculate 
the  average  velocity  for  each  calibration  point. 

15.  On  the  right  side  of  the  screen,  enter  the  number  of  calibration  points  to  be 
taken  (usually  8). 

16.  For  each  calibration  point,  enter  dP  and  the  press  enter  or  click  mouse 
outside  of  dP  box.  The  IFA  software  should  have  calculated  a  velocity 
from  the  dP  reading.  Press  [Acquire]  to  record  calibration  point.  Record 
data  on  calibration  sheet. 

17.  After  ah  calibration  points  have  been  taken,  close  graph  and  press  [Next 
Screen]. 
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18.  Select  [Curves]. 

19.  At  bottom  of  screen,  in  the  “fit”  block,  select  King’s  Law.  Record  King’s 
Law  coefficients  on  calibration  sheet. 

20.  Open  new  calibration  file  in  Notepad.  Change  probe  number  on  lines  2 
and  6  to  appropriate  number  and  save. 
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APPENDIX  C.  HOT-FILM  MEASUREMENTS 


1 .  Create  new  folder  with  today’s  date. 

2.  Copy  probe  calibration  files  into  new  folder. 

3.  On  main  menu,  select  “Acquisition”  then  “Probe  Table.” 

4.  In  Experiment  box,  select  “Get  File.”  Open  newly  created  folder.  Name 
new  file,  e.g.  90PE.R0001  for  90%  speed  peak  efficiency. 

5.  Click  on  “Add  Probe.” 

6.  Select  calibration  file  for  channel  1  transducer. 

7.  Repeat  steps  5  and  6  for  channel  2. 

8.  Ensure  both  channels  have  asterisks.  If  not,  select  channel  and  press  “Tag 
A/D.” 

9.  Click  “Next  Screen.” 

10.  Update  “Atm  Pressure”  and  various  units. 

11.  Set  sample  rate  to  100,000  Hertz,  size  to  8  kpts/ch. 

12.  Click  “Next  Screen.” 

13.  Once  TCR  speed  steady,  click  “Trigger”  to  record  data. 
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APPENDIX  D.  MATLABCODE 


A.  HW.M 

%  This  script  will  calculate  the  average  velocity  and 
%  turbulence  intensity  of  up  to  four  probes  by  first 
%  sending  raw  data  into  matrix  "a"  and  then  expanding 
%  the  matrix  with  the  calculated  instantaneous  velocities. 

%  The  velocities  from  probes  1,2, 3, 4  will  be  placed  in 
%  in  columns  6, 7, 8, 9  respectively,  irregardless  of  the  number 
%  of  probes  used. 

%file=input('Enter  probe  data  file  with  single  quotes  in  this  format:  filename. Exxx') 

%  format  for  file  is  T3JUL05.E0005' 

%  kings  law  coeff  must  be  in  a  matrix  with  probes  as  rows,  (a  b  n) 

%  for  example  k=[.87726  .31691  .5] 

king=[-23. 44088  20.06122  0.11]; 

a=dlmread('90  0PEN.E0003'); 

for  n=l:length(a) 
a(n,2); 

expo=l/king(l,3); 

E2=a(n,2)*a(n,2); 

a(n,3)=((E2-king(  1 , 1  ))/(king(  1 ,2)))Aexpo; 
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end 


%  Calculate  the  mean  velocity 

avevel_l=mean(a(:,3)) 

sd_l=std(a(:,3)); 

%  And  the  turbulence  intensity  from  the  standard  deviation 

turb  int_  1  =(  sd_  1  *  1 00  )/ave  ve  1_  1 

%  FFT 

time=a(:,l); 

lengtht=length(time); 

deltat=(time(lengtht,  1  )-a(  1 , 1  ))/lengtht 

plot(time,a(:,3)) 

Y=fft(a(:,3)); 

Y(l)=[]; 

nn=length(Y); 

power=abs(Y(l:nn/2)).A2; 

nyquist=l/2; 

%  freq=(l:nn/2)/(nn/2)*nyquist; 
freq  =  (l:nn/2)/(nn*deltat); 
figure 

loglog(freq, power) 
xlabel('frequency') 
ylabel('Power') 
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B.  FFTAUTO.M 

%  Script  to  plot  the  FFT  and  Autocorrelation  of  a  hotwire  signal 

close  all 
clear  all 

a=dlmread('90_PE.E0005'); 

%time=a(l:  1024,1); 

%voltage=a(  1 : 1 024,2); 

time=a(:,l); 

voltage=a(:,2); 

figure)  1) 

plot(time, voltage) 
xlabel('time') 
ylabel('voltage') 
title('Raw  Signal') 

N=length(voltage)- 1 ; 
deltat=.  00001; 

Y  =  fft(voltage); 

Y(i)=[]; 

n=length(Y); 

power  =  abs(Y(l:n/2)).A2; 
nyquist  =  .5; 
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%freq  =  (l:n/2)/(n/2)*nyquist; 
freq  =  (l:N/2)/(N*deltat); 
figure(2) 

loglog(freq,power) 

%semilogy  ( freq, power) 
xlabel('freq') 
ylabel('power') 
title('Periodogram') 

figure(3) 

ave_volt=mean(  voltage); 
volt=voltage-ave_volt; 

%Plot  the  fluctuating  voltage 

plot(time,volt) 

figure(4) 

%No.  of  points  to  take  the  autocorrelation  over 
m=8192; 

[C,lags]=xcorr(volt,m,'coeff); 

%plot(time,C) 

plot(lags,C) 

figure(5) 

auto_time=deltat*lags; 
plot(auto_time,C) 
grid  on 
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APPENDIX  E.  95%  SPEED  PEAK  EFFICIENCY  AND  NEAR  STALL 
AUTOCORRELATION  AND  CONTOUR  PLOTS 


Figure  47.  Auto-Correlation  of  Hot-Film  Data  at  90%  Speed  Near  Stall 


At  «  0.002 
Ti  =  0.019 
Uavg  =  65.8 

u’2  =  (Ti  *  Uavg)2  =  1.563 
k  =  1.5*(u’2)  =  2.3445 
L  =  At  *  Uavg  =  0.1316 
s  =  (C0  75  *  kL5)/(K  *  L)  =  1 1.206 


(From  above  graph) 

(Turbulence  intensity  from  hwl  .m) 
(Average  velocity  from  hwl.m) 

(Turbulence  length) 

(C  =  0.09,  k  =  0.4) 
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Figure  48.  Kulite  1  at  95%  Speed  Peak  Efficiency 
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Figure  49.  Impact  Probe  at  95%  Speed  Peak  Efficiency 

58 


i - 1 — - 1 - — - 1 - 1 - 1 - = — i  - 1 - 1 - r 


0.02  0.04  0.06  0.08  0.1  0.12  0.14  0.16  0.18  0.2  0.22 


Time 


log(Frequency)  log(Frequency) 


4.5 

4 

3.5 

3 

2.5 

2 

Figure  50.  FIot-Film  1  at  95%  Speed  Peak  Efficiency 
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Kulite  1  at  95%  Speed  Near  Stall 
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Figure  51. 
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Figure  52.  Impact  Probe  at  95%  Speed  Near  Stall 
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Figure  53.  Flotfilm  1  at  95%  Speed  Near  Stall 
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APPENDIX  F.  KULITE  CONTOUR  PLOTS  FOR  90%  AND  95% 
SPEED:  STEAM  INDUCED  STALL 


Figure  54.  Kulite  2  at  90%  Speed  Steam-Induced  Stall 


Time 


Figure  55.  Kulite  3  at  90%  Speed  Steam-Induced  Stall 
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Figure  56.  Kulite  4  at  90%  Speed  Steam-Induced  Stall 
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Figure  58.  Kulite  6  at  90%  Steam- Induced  Stall 
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APPENDIX  G.  CFX  RESULTS 


CFX  at  90% 

Run  # 

nidot  ave 

Pt  in  abs 

Pt  out  abs 

Pi 

Tin 

Tout 

Tao 

Efficiency 

TCR  0 

7.2700 

101329 

128781 

1.2709 

300.0 

324.6 

1.0822 

86.24 

TCR  2 

7.2518 

101329 

130598 

1.2888 

299.9 

325.6 

1.0857 

87.72 

TCR  4 

7.2259 

101329 

132433 

1.3070 

299.9 

326.7 

1.0894 

88.93 

TCR  6 

7.1876 

101329 

134278 

1.3252 

299.8 

327.7 

1.0932 

89.85 

TCR  8 

7.1318 

101329 

136118 

1.3433 

299.7 

328.8 

1.0972 

90.51 

TCR  10 

7.0529 

101329 

137926 

1.3612 

299.5 

329.8 

1.1012 

90.98 

TCR  12 

6.9486 

101329 

139714 

1.3788 

299.3 

330.8 

1.1053 

91.29 

TCR  14 

6.8140 

101329 

141520 

1.3966 

299.0 

331.8 

1.1096 

91.35 

TCR  16 

6.6255 

101329 

143403 

1.4152 

298.6 

333.1 

1.1153 

90.51 

TCR  17 

6.4686 

101329 

144474 

1.4258 

298.3 

334.0 

1.1196 

89.19 

TCR  18 

3.7401 

101313 

136621 

1.3485 

295.4 

330.3 

1.1179 

75.62 

CFX  at  95% 

Run  # 

nidot  ave 

Pt  in  abs 

Pt  out  abs 

Pi 

Tin 

Tout 

Tao 

Efficiency 

TCR  0 

7.5392 

101329 

131541 

1.2981 

300.6 

328.1 

1.0917 

84.39 

TCR  2 

7.5324 

101329 

133391 

1.3164 

300.6 

329.1 

1.0950 

86.01 

TCR  4 

7.5220 

101329 

135251 

1.3348 

300.5 

330.0 

1.0982 

87.61 

TCR  6 

7.5093 

101329 

137189 

1.3539 

300.5 

331.1 

1.1020 

88.69 

TCR  8 

7.4895 

101329 

139124 

1.3730 

300.5 

332.2 

1.1056 

89.73 

TCR  10 

7.4604 

101329 

141084 

1.3923 

300.4 

333.3 

1.1096 

90.52 

TCR  12 

7.4166 

101329 

143064 

1.4119 

300.3 

334.5 

1.1138 

91.02 

TCR  14 

7.3481 

101329 

145019 

1.4312 

300.1 

335.6 

1.1181 

91.29 

TCR  16 

7.2485 

101329 

146922 

1.4499 

299.9 

336.7 

1.1225 

91.39 

TCR  17 

7.1863 

101329 

147872 

1.4593 

299.8 

337.2 

1.1248 

91.35 

TCR  18 

7.1141 

101329 

148838 

1.4689 

299.6 

337.8 

1.1273 

91.20 

TCR  20 

6.8936 

101329 

150955 

1.4897 

299.1 

339.3 

1.1344 

89.78 

TCR  22 

3.3971 

100634 

140799 

1.3991 

294.7 

336.8 

1.1429 

70.48 
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